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Quasi-one-dimensional quantum structures with spectra scaling faster than the square of 
the eigenmode number (superscaling) can generate intrinsic, off-resonant optical nonlin¬ 
earities near the fundamental physical limits, independent of the details of the potential 
energy along the structure. The scaling of spectra is determined by the topology of the 
structure, while the magnitudes of the transition moments are set by the geometry of the 
structure. This paper presents a comprehensive study of the geometrical optimization 
of superscaling quasi-one-dimensional structures and provides heuristics for designing 
molecules to maximize intrinsic response. A main result is that designers of conjugated 
structures should attach short side groups at least a third of the way along the bridge, 
not near its end as is conventionally done. A second result is that once a side group is 
properly placed, additional side groups do not further enhance the response. 

Keywords: nonlinear response; hyperpolarizability; fundamental limit; quantum-confined 
systems; quantum wires; sum rules 


1. Introduction 


liquidaiZI’EfniilcJ, liquid cryst 
mesoscopic solid state wire; 



Nonlinear optics is the study of quantum systems with polarizabilities that are 
nonlinear functions of external elecLmmagnetic fields. Nonlinear optical phenom¬ 
ena include harmonic generation Jj’BtJ, optical phase conjugationclfl and four- 
wa^e mixingpEl, optical bistabihtyEm, ultrafast optics 119 , and waveguide 
ii i 2* A/iof^oio w ith suitable polaiizabilities have been sought in solidal3 

and gasesE3’E2l, as well as in ohotonic crystals, 
and artificial systemsEZl’Ej’OLy. Applications 
center on electric-field control of light for ultrafast optical signal modulation using 
the third-rank optical susceptibility tensor Xijk and all-optical, light by light con¬ 
trol for high-speed fiber-optic switching using the fourth rank optical susceptibility 
tensor Xijki • Ultrafast effects, which rely on virtual transitions among electronic 
states, arise from the far-from-resonance response to an electric field. A prime ob¬ 
jective of research in nonlinear optics is the discovery of the fundamental origins of 
large intrinsic responses that may be realized in macroscopic materials and artificial 
structures. 


l 
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Fund am en talg uant 11 m mechanics limits the intrinsic nonlinearities of any molec¬ 
ular structureLllSH’O. Until 2006, a factor of thirty gap existed between theoret¬ 
ical and experimental values of the first hyperpolarizabihtyf!3. The existence of the 
gap suggested that conventional molecular design rules required modification to 
enhance the intrinsic nonlinear opticalresponse. Advances in the design and syn¬ 
thesis of conjugated chain polyynesHOES, and extended conjugation and substi¬ 
tuted donor-acceptor cyanoethynylethene moleculesBUj produced molecules with 
larger third-order nonlinear optical susceptibilities, as did modulated conjugation 
molecular structures 0000 

For off-resonant phenomena, the eigenstates and spectra of the system are the 
critical parameters that control the hyperpolarizabilities. The first comprehensive 
study of the relation between spectral properties and the nonlinear optical response 
was a Monte Carlo calculations which showed that a necessary condition for a large 
intrinsic response is that the energy spectrum should scale faster than the square 
of the eigenmode number, so-called superscaling spectras. Superscaling spectra 
resemble particle-in-the-box energies, with periodic boundaries separating randomly 
positioned wavenumbers. It has been shown that the potential in donor-acceptor 
configurations can produce such spectra. 

Donor-acceptor structures with superscaling spectra require another design rule 
to achieve maximum response, one that specifies the relative distributions of the 
lowest eigenstates lying above the Fermi level in a multielectron system. Near the 
quantum limits, the response for the first hyperpolarizability is dominated _hy the 
contributions from only three eigenstates, the so-called three-state conjecturel!3. For 
the second hyperpolarizability, the response is dominated by four eigenstates. The 
optimum transition moments are those creating spatial charge separation among 
the ground and excited states, leading to a large change in the dipole moments 
x nn = x nn — xoo for n = 1,2, and with an additional requirement that the ge¬ 
ometry maximize the alignment of the lowest order off-diagonal transition moment 
contributions £ 015 ^ 02,^12 to the response. 

This paper is the third in a series of studies of the optimum topology and ge¬ 
ometry for quasi-one-dimensional structures, such as nanowires or donor-acceptor 
molecules, by analyzing the eigenstates and spectra of one electron quantum graph 
models with bare and dressed edges. The first study showed that dressed quantum 
graphs could achieve nonlinearities near the fundamental limited, while the sec¬ 
ond revealed the optimum topologiesej for producing such spectra. In this paper, 
we show that the optimum geometries for systems with side groups and optimum 
spectra are those that generate the requisite change in dipole moments while main¬ 
taining a large spatial overlap among the transition moments between the lowest 
three or four eigenstates. The quantum graph models suggest that a general rule 
for designing real structures is to place the side groups at specific locations along 
the chains, but far from the ends. Such configurations create eigenfunctions pro¬ 
ducing large changes in the dipole moments in the presence of an external optical 
field while at the same time allowing good spatial overlap for large off-diagonal 
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transition moments. 

Quantum graphs are quasi-ID systems to which electron dynamics are con¬ 
fined along the edge of a network structure which can be thought of as a 
branched nano-wire structure, or a quasi-linear molecule, such as a donor-acceptor, 


with skk 
models@ 
tical waveguides 
graphene, and carbon m 
models of quantum chao: 



Quantum graphs were first studied as tractable molecular 


systemsO, op- 
and fullerines, 


,and have been i: 
quantum win 



d as models of mesoscor 
excitations in fractals t 


il. Quantum graphs are also exactly solvable 
rave recently been success- 


Microwave network 

fully used to experimentally simulate quantum graphsEfk 

Quantum graphs were recently introduced into nonlinear opticsLZIH’L^LyO’O as 
model systems exhibiting the optimum spectra for large intrinsic nonlinear optical 
responses, for specific geometries. In this paper, we exploit quasi-one-dimensional 
quantum graphs in a comprehensive study of the optimum geometries of nonlinear 
optical structures with particular focus on the placement of a defect or side group 
along a main chain such that the spatial distribution of the lowest eigenfunctions 
produces transition moments that optimize the response. 

All of the graphs analyzed in this paper are of the form of a wire supporting 
a S function potential simulating a defect, a side prong simulating a side group, 
and combinations of these. Typical donor-acceptor molecules contain two identical 
groups on each end, one nearly collinear with the main chain and the other at 
an angle to it. Dressed and prong wire graphs simulate the geometry of donor- 
acceptor molecules, but their dynamics are obviously an idealization of a many 
electron molecular system. However, the single electron model shows explicitly the 
form and shape of the states and spectra that real systems must possess. Since 
chemists are expert at shaping states and spectra with strategic placement of side 
groups along main chains, the heuristics described here should prove quite useful to 
them in their pursuit of molecular designs with the geometry and topology required 
to bridge the gap. 

This paper is organized as follows. Section || briefly reviews the calculation of the 
first and second hyperpolarizabilities of quantum graphs with arbitrary shape and 
topology. Section [| discusses graphs dressed with S functions. Section || discusses 
the graphs with prongs only. Section [| discusses dressed graphs with prongs. We 
show that S graphs with a positive potential are spectrally similar to prong graphs 
for the lowest energy eigenfunctions on the graphs. We derive a design heuristic 
for the graphs to generate the optimum transition moments. Section presents our 
conclusions and discusses the applications of the heuristics to molecular design. A 
set of appendices reviews two other kinds of graphs, those with square wells or steps 


in Appendix A, and those with slant wells in Appendix B 
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2. Nonlinear optics of planar quantum graphs 

The calculation of optical nonlinearities for quantum graphs is articulated ir 
In brief, the eigenfunctions and spectra are determined by solving the nonrelativis- 
tic wave equation for the edges of the graph and applying conservation conditions 
at the vertices. The spectra and transition moments are then used in a sum over 
states to compute intrinsic second and third order nonlinearities Ed. Tensor proper¬ 
ties are easily studied usin^ spherical representations of the Cartesian tensors /3ijk 
and 'jijki , as described in L3. All quantum graphs have superscaling spectra, with 
wavenumbers contained within root boundaries spaced more or less equally. The 
following subsections summary our standard method. 


2.1. A quasi-ID dynamical system 

The dynamics of an electron on a quantum graph are described by a self-adjoint 
Hamiltonian operating on the edges of the graph, with complex amplitude and 
probability conservation (hereafter referred to as flux conservation throughout the 
paper) at all internal vertices and fixed, infinite potentials at the termination ver¬ 
tices (where the amplitude vanishes). The physics of the eigenfunctions and their 
spectra have been previously described, along with a suitable lexicography for de¬ 
scribing the union operation for creating eigenfunctions from the edge functions 
that solve the equations of motion for the HamiltonianEd. 

The graph is specified by the location of its vertices and the edges connecting 
the vertices. Figure [j] details the notation and configuration of a graph. A set of 
vertices with arbitrary locations in the 2D plane but fixed connections specifies a 
topological class of graphs. For a fixed topology, the variation of vertex locations 
specifies various geometries for the graph. Since motion is confined to the graph 
edges and is continuous at each vertex, the energy spectrum depends only on the 
edge lengths and the boundary conditions, ie, the topology. Spectra are quasi¬ 
quadratic in state number, ie, superscaling. This is the first requirement for a large 
response from a nonlinear optical structure. 

The edge lengths and angular positions determine the projections of electron 
motion onto a fixed, external reference axis. The projections summed over all edges 
yield the transition moments required to compute the tensor elements of the hy¬ 
perpolarizability tensors. Regardless of how the axes used to define the vertices are 
chosen, the various tensor components may be used to assemble any component in 
a different frame by using the rotation properties of the tensors. 

The study of the nonlinear optical properties of a specific graph topology re¬ 
quires solving the graph for its eigenfunctions and spectra as functions of its edge 
lengths and using them to compute a set of transition moments for the graph from 
which the hyperpolarizability tensors may be computed. 
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Fig. 1. A four-edge quantum graph. Each edge has its own longitudinal coordinate Si ranging 
from zero to L^. The projection x(si) of an edge onto the x — axis is measured from the origin 
of the coordinate system attached to the graph (and not to the beginning of the edge itself). For 
example, x(s i) = s i cos6h while x(s 2 ) = L\ cos6h + S 2 cos 62 and so on. 


2.2. States and spectra 

A general quantum graph is solved for its spectra and states once the Hamiltonian 
has been specified. Then, the Hamiltonian is used to compute a set of edge func¬ 
tions, </>^( Si ) that are solutions to the Schrodinger equation with the same set of 
eigenvalues E n on each edge, but with specific boundary conditions determined by 
their connections at the vertices. Next, the edge functions are used to construct 
eigenfunctions of the Hamiltonian for the graph through union process that reflects 
a direct sum Hilbert space over the edges lJ 

V’n(s) = U f =1 4> l n (Si) (2.1) 

The union is defined such that the eigenfunction is continuous at every vertex of the 
graph, while the probability current is conserved at each vertex and thus through¬ 
out the graph. These two boundary conditions guarantee that the eigenfunctions 
are complete. They also generate the relationships among the edge amplitudes re¬ 
quired to compute the spectrum of the graph. Solutions of the amplitude equations 
resulting from the boundary conditions exist only if the determinant of the matrix 
of coefficients vanishes. This condition produces the secular or characteristic equa¬ 
tion for the graph and determines the eigenvalues k n and the exact energy spectrum 
E n . Since the boundary conditions in the elementary QG model are independent 
of the angles the edges make with respect to one another, the secular equation is 
independent of angles and depends only on dimensionless parameters k n ai. For a 
given configuration of vertices, the distance between them and the rules by which 
they are connected, i.e., the topology of the graph, determines the energy spec- 
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trum. Topologically different graphs with identical geometries have different energy 
spectra. In this way, the graph topology has a large impact on the nonlinear optical 
response. 

Except for bent wires and closed loops, the secular equation of a graph is gen¬ 
erally a transcendental equation. Accurate solutions are easily found numerically. 
From these, the relative internal amplitudes may be calculated. Normalizing the 
eigenfunction produces the states required to compute the transition moments. 

The transition moments are sums (not unions) over edges of the following form: 

E rai 

Xnm = E / ( s *)0m( s *) x i s i) ds i ( 2 - 2 ) 

i= 1 Jo 

where <fi l m (si ) is the normalized wave function on the i th edge that obeys the bound¬ 
ary conditions for the graph and x(si) is the x-component of s*, measured from the 
origin of the graph (and not of the edge), and is a function of the prior edge lengths 
and angles. 

The transverse wavefunction is not calculated and not essentiaLin this model. 
The transverse states do not contribute to the hyperpolarizabilitiesO. But there are 
residual effects from the transverse state in the sum rules, as has been previously 
discussedOO. 

The process for calculating the spectra and transition moments of elementary 
QG’s can be summarized as follows: (1) select a particular graph topology, speci¬ 
fying the number of vertices and the connecting edges, (2) generate a random set 
of vertices, and calculate the lengths of the edges and the angles each makes with 
the x-axis of the graph’s coordinate system, (3) solve the Schrodinger equation on 
each edge of the graph, and (4) match boundary conditions at the vertices and 
terminal points. This results in a set of equations for the amplitudes of the wave- 
functions on each edge. The solvability of this set requires that the determinant of 
the amplitude coefficients vanishes, leading to a secular equation for the energies. 
The transition moments x nm and energies may be used to compute the first and 
second hyperpolarizabilities of any graph specified by a set of vertices, as described 
next. 


2.3. Hyperpolarizabilities of planar quantum graphs 


For the quasi-ID problem the tensors are indexed in the (x,y) directions. The full 
tensor expressions are given as a sum over states. The first intrinsic hyperpolariz- 
ability tensor for 2D graphs may then be written as 


fiijk — 


P 


Pn 


I) £ 


3/4 


'CoJlmtimO 


e n &r 


(2.3) 


where and e n are normalized transition moments and energies, defined by 


T = 

Snm 


E, 


„max ' 

01 


^77, - 


n0 


Ei o 


(2.4) 
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with A 1 ^ = x and A 1 2 ) = y , E nm = E n — E m , and where 


f h 2 \ 1 ^ 2 

max _ ( lv \ 

01 ~~ \2mE 10 ) ' 


(2.5) 


r 01 


represents the largest possible transition moment value of roi0 According to 


Eq. (E4), eo = 0 and e\ = 1. fiijk has been normalized to its maximum value, is 
scale-invariant, and can be used to compare molecules of different shapes and sizes. 
Similarly, the second intrinsic hyperpolarizability is given by 


fijkl — 


E 

i n.m.l 


'Zon&mZLti 


Cn^m^l 


-E 


ti £3 £J tk 
SOnSnOSOmSmO 


eter, 


( 2 . 6 ) 


The maximum valuesmsed for normalization are derived from the Thomas-Reiche- 
Kuhn sum rulesLd’HHo, and are 


and 


Pmax = 3 1/4 


eh 

nM 2 


TV 3 / 2 


E 


,7/2 

10 


7ma 


= i( e X\ 


N 2 


\m 2 J Ef { 


10 


(2.7) 


( 2 . 8 ) 


where E \o is the ground state to first excited state energy gap, and N is the number 
of electrons. The second hyperpolarizability normalized this way has a largest nega¬ 
tive value equal to —(1/4) of the maximum value. Hyperpolarizabilities normalized 
this way enable direct comparisons of the intrinsic response without regard to size. 

The fundamental limits for fallshort of unity for every system calculated 
from a potential energy functionOcj’O and appear to asymptote to a value of 
0.7089. For 7 ^^, the apparent maximum and minimum limits are sixty percent 
of the values predicted by Eq. fl2.8| ). For the rest of this paper, we use the term 
fundamental limits to mean the apparent fundamental limits. 


3. Dressed graphs 

Figure [| shows a bare wire graph and three variants, each with a modification of the 
bare wire to introduce finite 5 potentials. The addition of one or more 5 functions 
to a wire creates a source or sink of probability fluxEl. This has a profound effect on 
the spectrum at the lower end, and on the shape of the ground state. Both effects 
are expected to cause an enhancement in the hyperpolarizabilities. The details of 
how this comes about are described in this section. 
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Fig. 2. The bare wire graph and variants with one, two, and three S potentials, respectively. The 
probability flux at the position of each 5 potential is discontinuous by an amount determined by 
the strength of the potential. 


3.1. Bare wire 

Figure [|, upper left panel, shows a bare wire quantum graph consisting of an edge 
supporting electron dynamics, bounded at either end by an infinite potential. The 
electron behaves as a particle-in-a-box, with /3 XXX = 0, and ^y xxxx = —0.126. How¬ 
ever, the graph is a quasi-ID object, and if the wire is bent or is comprised of 
two edges with a nonzero angle between them, (3 XXX is nonzero, and 7 xxxx changes 
slightly. The first seven eigenfunctions for this graph are shown in the left panel 
of Figure ||. For the bare wire (one edge), the centrosymmetry keeps /3 XXX at zero. 
The probability flux in the well varies from one edge to the other in a continuous, 
conserved fashion. The probability that the electron is on the edge is always unity. 

We now examine variations of this graph that have been modified to include 
finite S potentials. The modifications create interesting quantum graph structures 
that are straightforward to solve, have more complicated energy spectra, and have 
eigenfunctions whose shapes are altered in such a way as to create giant enhance¬ 
ments of the optical nonlinearities by localizing the ground state to the 5 potentials, 
while pushing the first few excited states out from the ground state. 


3.2. Wire with one 6 potential 

The upper right panel in Figure [ 2 ] shows a wire to which a finite potential of the 
form V(x ) = ( g/L)S(x — a) has been placed at location x — a. The nonlinear optics 
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of this graph have been discussed extensively^. The S function is a source or sink 
of probability flux, creating a slope discontinuity in the eigenfunctions of the wire. 
The right panel in Figure || shows the localization of the ground state wavefunction 
to the position of the potential, with little to no change in the localization of the 
first few excited states. This delocalization of the ground state from the first few 
excited states drives the nonlinearities from zero to values near the fundamental 
limits of the off-resonance hyperpolarizabilitiesLT This is even more evident in for 
Figure | for ~ fxxxx . 
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Fig. 3. Lowest eigenfunctions for a bare wire graph with zero potential (left) and one with a 
strong, negative 5 potential (right) positioned to maximize (3 XXX . 
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Fig. 4. Lowest eigenfunctions for a bare wire graph with zero potential (left) and one with a 
strong, negative 5 potential (right) positioned to maximize 7 xxxx . 


Figure shows the hyperpolarizabilities varying with position of the S function 
on the wire at different strengths g. The spectrum changes when the strength eps 
smaller than a critical strength g c determined by the location of the S functions. 
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For g < g c < 0, the ground state is negative, and the deep well pulls the lowest 
eigenfunctions to the S function, creating a large overlap and generating a large 
response for /3 XXX , even as the S function moves toward one end of the wire. For 
positive g, /3 XXX gradually goes toward zero as the S function moves toward either 
end. In both cases, /3 XXX vanishes when the S function is at either end, due to the 
centrosymmetry of the resultant wire. For j xxxx , large values are obtained as the S 
function moves toward the middle of the wire when the well is very deep. 




0.0 0.2 0.4 0.6 0.8 1.0 

Position of 5 function 


Fig. 5. Variation of (3 X xx and r y X xxx with position along the wire of the S function for various 
strengths. 


Figure |] shows the eigenfunctions for a graph with a barrier and one with a well 
of equal potential. The former looks like a bare wire, with near-zero response, while 
the latter achieves optimum response near the fundamental limits. To the extent 
that a negative S function mimics a quantum defect causing disruption of electron 
flow, the Figure suggests that the presence of a simple defect on a linear structure 
is sufficient to displace the lowest eigenfunctions to create a large nonlinear optical 
response. 

Figure [ 7 ] shows the relationship between eigenfunction localization and large 
response. In the Figure, the electron density across the wire is shown as a shaded set 
of bars corresponding to the first five eigenfunctions, positioned at three locations 
of the S function on a half-wire. Highest density has the greatest shading. Also 
shown is the value of /3 XXX as a function of the position of the S function. The 
lowest response occurs when the S function (a well, in this case, since g = —3.73) is 
near one end of the wire or at the middle of the wire. The greatest shift in charge 
density between the ground and lowest excited states occurs when the 5 function 
is located a little over a quarter of the way along the wire. The electron density is 
delocalized, and /3 XXX is small. 

The analysis of the one S graphs suggests that quantum structures with defects 
located at precise points along a wire can produce large enhancements of the non¬ 
linear optical response. A positive S function yields a large response if it is about 
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Fig. 6. First seven eigenfunctions for the (5 graph with positive strength g, near one edge (left), 
with a small /3 XXX , and equal but negative strength g, with a large /3 XXX . 


0.0 0.2 0.4 0.6 0.8 1.0 



Distance x along the wire 


Fig. 7. Electron charge density for the first five eigenfunctions (with the ground state at the 
bottom of each band) along the length of the 15 graph when the S function with g = —3.73 is 
placed at three different locations. 


midway between the center and edge of the structure; a negative S function has its 
greatest effect closer to one end of the structure. Modulation of the location of the 
defect by an external control electric or magnetic field would cause a huge change 
in the nonlinear optical refractive index, as is required for an optical modulator or 
bistable optical switch. 
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3.3. Two ds 

A single S function on a wire causes a large change in the nonlinear optical response, 
generating intrinsic nonlinear it ies near the apparent fundamental limits. The lower 
right panel in Figure shows a wire with two S functions on it. This graph was 
previously solveden. The eigenfunctions for the graphs with the largest /3 XXX and 
Ixxxx are shown in Figure [| 



Fig. 8. First seven eigenfunctions for the two 5 graph with the largest (3 X xx and jxxxx- 


Figure [ 9 ] displays the dependence of the hyperpolarizabilities when one S func¬ 
tion is fixed and the other scans across the wire. Figure 10 displays the dependence 
of the hyperpolarizabilities when one S function is fixed near the left edge (a = 0.14) 
with fixed strength g 1 = 5 and the other is positioned at a fixed distance b from 
the right edge, but for several values of # 2 - These are essentially vertical slices in 
the prior Figure and show that the combination of a barrier near one end and a 
well at the other results in a large response for either hyperpolarizability, a toy 
donor-acceptor model for a quantum wire molecule. 

show the dependence of the hyperpolarizabilities on the strengths of 


Figure 


11 


the two 5 functions. 

The results of this work indicate that an additional defect has little effect on the 
nonlinear optical response, once a single defect is in place at the proper location. 


3.4. Three ds 

The lower right panel of Figure || shows three S functions on a wire. This graph was 
previously studied in detailed. The eigenfunctions for the graphs with the largest 
/% xx and ^f xxxx are shown in Fig [l|. 

Prior analysis showed that the hyperpolarizabilities approached their maximum 
values with one S function with the proper placement and strengthen . The max¬ 
imum values were just shy of the theoretical maxima implied by the theory of 
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Position of 2nd § function 



Position of 2nd 5 function 


Fig. 9. Line plots for f3 xxx at a — 0.14, g\ — 5. This corresponds to a 8 barrier near the left end 
of the wire. The various curves show how (3 XXX and r y xxxx change as the potential from the second 
8 function is moved toward the other end and at various strengths corresponding to both barriers 
and wells. 




g 2 g 2 

Fig. 10. Line plots for / 3 XXX for a = 0.14, g\ = 5 as the second (5 function is fixed at several 
distances from the right edge, with varying strengths. 


fundamental limits. The addition of a second S function produced a very slight 
increase in both hyperpolarizabilities, and we see that the three S wires have hyper- 
polarizabilities that are essentially at the fundamental limits. We conclude that it 
is not necessary to interrupt the flux on a donor-acceptor chain with more than one 
S function or prong. This is in agreement with prior analysis showing that only a 
few parameters are required to achieve near-maximum nonlinear it ies, despite that 
this poorly defines the potential 7 . 

Again, the presence of an additional defect does little to enhance the nonlinear- 
ities, but it also does not negatively affect them. This would suggest that a wire 
with many randomly placed S defects, would likely have at least one in a position 
generating a large response, and that the presence of the rest of the random defects 
would have minimal positive or negative impact. 
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Fig. 11. Contour plot for the two S graph of (3 XXX and r ) xxxx for fixed S locations but variable 
strengths. 



Fig. 12. First seven eigenfunctions for the three 5 graph with the largest /3 XXX and 7 xxxx 


4. Quantum structures with side groups 

Figure |l^ shows a bare wire graph and three variants, this time with prongs at¬ 
tached. The addition of a prong at a junction to a bare wire creates a 3-star graph, 
previously studied in detail for its quantum mechanical properties!^ and solved for 
its nonlinear optical propertiesllj. One of the variants has the geometry of the series 
of donor-acceptor structures^! The presence of a prong creates a pathway for charge 
flow away from the main chain. This changes the spectrum to a quasi-random set of 
energies lying within fixed root boundaries equal to the energies of the bare wire. It 
also changes the shape of the eigenfunctions, affecting the change in dipole moment 
and the off-diagonal transition moments. The proper choice of a prong length and 
position can lead to a large enhancement in the hyperpolarizabilities. 
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L 



c 



Fig. 13. The bare wire graph and a variant with one prong, one with two prongs, and another 
with two prongs but with a geometry suggestive of a donor-acceptor molecule. The probability 
flux at the position of the prongs is conserved but the prong carries away some of the probability 
flux, creating an apparent discontinuity across the straight part of the wire at the prong location. 


4.1. Wire with one prong 


The upper right panel of Figure 13 shows a bare wire graph to which a single prong 
has been attached instead of a S function. This graph has a star topology. Figure 
0 show the variation of the hyperpolarizabilities with position of the prong on the 
wire at a number of prong lengths. The behavior resembles that of the S graph when 
g > 0. There is no analog to the negative g S graph for the prong graph because 
the energy spectrum of star graphs is always positive. 

Figures [l5| compares the first seven eigenfunctions of the bare wire graph (left) 
with those of the one prong graph with the largest /3 XXX . The enhancement induced 
by the existence of a short prong arises from the change induced in the distribution 
of electron flux in the ground state, compared to the relatively unchanged shape of 
the first few excited states. 

Figure 0 is a 3D plot of the edge vectors for the 1 prong graph for the first 
three eigenfunctions. The two edges along the x-axis are of equal length, producing 
a degeneracy in the first excited state, as shown in the Figure. One degenerate state 
has a nonzero amplitude on all edges, while the other has exactly zero flux on the 
prong. The nonlinear it ies for this graph are well below the maximum values for this 
topology. 

Figure 0 is a 3D plot of the edge vectors for the 1 prong graph for the first 
three eigenfunctions when the two edges along the x-axis are irrationally-related, 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Position of prong on wire Position of prong on wire 


Fig. 14. Variation of /3 X xx and r ) xxxx with position of the prong for various prong lengths. 



Fig. 15. Lowest eigenfunctions for a bare wire graph (left) and one with a prong of length b = 0.08 
positioned to maximize /3 X xx This configuration has a large r ) xxxx as well. If the prong length is 
decreased to b = 0.05, then 7 xxxx —> 0.300, its maximum for the one prong graph, while (3 XXX 
drops slightly. 


producing a nondegenerate spectrum. The fi xxx for this graph is at the maximum 
values for this topology. The effect of the prong is to cause a large slope discontinuity 
along the edge, localizing the lowest order states and leading to a large response. 

How similar are prong graphs and S graphs with g > 0? The secular equation 
for the one prong graph may be written as 


cot ka + cot kc + cot kb = 0. 


(4.1) 


The secular equation for the one S graph may be written as 

, , 2 g 

cot ka + cot kc + —- = 0. 

kL 


(4.2) 


The equality of these two spectral equations would require an energy-dependent 
g(k) given by g{k) = ^ cot kb, but when k is near zero, g(k) and becomes 
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Fig. 16. Three-dimensional view of the edge states for the wire with one prong. In this graph, 
the two edges a and c are equal, and the prong is labeled b. The first excited state is two-fold 
degenerate, as shown. For this graph, (3 X xx is nonzero but small. The scale on the edges is arbitrary, 
as the intrinsic nonlinearities are scale-independent, but it is convenient to choose a + c = 1 along 
the wire axis. The vertical axis is labeled with the value of the amplitude of the normalized 
eigenfunctions. 


energy-independent. For b = 0.08, we find near-equality in the states and spectra 
when g ~ 6 . Thus prong graphs and S graphs are essentially equivalent in this case. 

To first order, the position for a prong on a wire that leads to maximum /3 XXX is 
evidently determined by maximizing the product \xoi\ 2 xu. Moving the prong from 
the edge toward the center pushes the ground and first excited states to opposite 
sides of the prong, which increases the change in dipole moment but decreases the 
ground to first excited state transition moment. The ideal position is somewhere 
in-between the edge and the center. The conclusion we can draw is that shaping the 
eigenfunctions in this manner will produce a large, intrinsic first hyperpolarizability. 
Similar reasoning holds for 7 xxxx when g > 0.. 

It appears that the genesis of large nonlinearities is a distortion of the eigen¬ 
functions along the main chain by the presence of the disturbance at the junction. 
The disturbance induces as phase shift in the part of the eigenfunction that lies 
along the main chain, producing a change in dipole moments sufficient to generate 
a larger /3 or 7 . To get a feel for how the one electron graph behaves as the prong 
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Fig. 17. Three-dimensional view of the edge states for the wire with one prong. In this graph, 
the two edges a and c are unequal, and all eigenfunctions are singlets. The prong interrupts the 
flux along the wire and generates the maximum (3 X xx for this topology. 


length is varied while its position is fixed, see the contour graph in Figure 18. It 



position of prong 


position of prong 


Fig. 18. Contour plot of \(3 XXX \ over the range of the position and length of the prong. 
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is evident that p is nearly zero for near-zero length prongs, regardless of where 
the prong is positioned. For a fixed prong position, say a = 0.37, increasing the 
prong length from zero causes a jump in p to its maximum for this topology, and 
continuing to increase the prong length beyond a value of about 0.08 causes p to 
drop again to a low value. For a fixed prong length, its ideal position is about a 
third of the way along the wire, while its least effective position is at either end, or 
at the center where the structure becomes centrosymmetric and f3 0. 

Moreover, the length of the prong attached to the main chain cannot be near¬ 
zero, nor can it get too large. To quantify this, abandon the stick model for a 
moment and consider the transition moments on the main wire of the chain and 
those on the side group. First, since the side group is effectively normal to the 
main chain, it doesn’t contribute to P xxx so we don’t want much of the transition 
moment to arise from this group. This means making it short so that an integral of 
the eigenfunction over that length is small compared to those on the main chain. 
That is precisely the reason we do not want a long side chain. But if the side chain 
is too short, the eigenfunction on the side chain develops a very large amplitude 
over a small distance, creating a steep slope at its junction with the main chain and 
causing the ground state on the main chain to completely localize to one side of the 
attachment point and the first few excited states to localize to the other side. This 
has the effect of driving the transition moments xoi and £02 to near-zero values, 
thus driving p and 7 to zero. 

Contributions to P (top) and the transition moments from which they arise 
(bottom) as the prong length is varied from small to large are shown in Figure |l9| . 
Note that it is Pn that contributes the most to P xxx . So we are most interested in 
how the contributing moments hi and xoi change with prong length b. 

Figure ( p0| ) shows exactly how the dependence in Figure ( |l9| ) arises. In the top 
left panel, the prong is very short and creates a large change in dipole moment, but 
the prong is so short that the ground and first excited state are localized to opposite 
sides of the prong, reducing their overlap transition moment xoi to near zero and 
effectively zeroing P xxx . As the prong length is increased (upper right panel), the 
change in dipole moment is reduced somewhat, but the overlap jumps, and the 
value of P xxx jumps to near its maximum for this topology. Continuing to increase 
the prong length keeps the overlap moments large but reduces the change in dipole 
moment to near zero, thus reducing p xxx again to a small value. The existence of an 
optimum prong length is highly suggestive of a design rule for side chain lengths on 
a donor-acceptor chain. This will be discussed shortly. The figure also shows how 
the three level model is nearly exact for P xxx for this graph. 


4.2. Two prongs 

The lower left panel of Figure [l3| shows a wire with two prongs. This graph is a 
back-to-back stamgraph, with two star vertices, but with the prongs perpendicular 
to the main wirecS. 
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Fig. 19. Contributions to (3 X xx (top) and the transition moments from which they arise (bottom) 
as the prong length is varied from small to large. 


Figure HI shows the eigenfunctions for the two prong graphs with the largest 
value of the hyperpolarizabilities. Figure displays the variations of f3 xxx and 
Ixxxx as each of the two prongs varies in size. One prong has been located near 
the center of the wire, while the other is midway between the center and the right 
edge. As is obvious from the eigenfunctions, short prongs are ideal for introducing 
a flux interruption along the wire, separating the ground state flux from those in 
the first few excited states. Figure [m] presents a parametric view of the variation 
with prong lengths. 
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Fig. 20. Edge states along the main chain for short prongs. 



Position x on wire 


Position x on wire 


Fig. 21. First seven eigenfunctions for the two prong graph with the largest (3 XXX and 7 xxxx 


4.3. Two prong wire with bent connector 

The lower right panel of Figure |l^ shows a two prong graph but with a bent wire 
connecting the two star nodes and with the geometry of the molecules shown in 
Figure [m[ These represent a subset of small cyanoethynylethene molecules in which 
donor-acceptor substitution was used to increase the second hyperpolarizability&S. 
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Fig. 22. Contour plot for the two prong graph of /3 X xx (left) and r ) xxxx (right) for fixed prong 
locations but variable lengths. 




0.00 0.05 0.10 0.15 0.20 0.25 

2nd prong length (c) 


Fig. 23. Variation of (3 XXX and r ) xxxx with the length of the second prong c (located at x = a + e) 
for various lengths of the first prong b (located at x = a), over a unit length wire. 


Figure g5| shows the variation in / 3 XXX and 7 xxxx as the link between the two 
bent wires (double bonds) is extended. If we assign values typical of single, double, 
and triple Carbon bonds, then the lengths of the fixed edges of the graph in the 
inset figure are calculable, as is the variable length of the linker as additional single- 
triple bond extensions are added. The result is a plot showing the existence of an 
ideal linker for (3 XXX (and a slightly longer one for 7 xxxx ). These results are again 
suggestive of a design rule for donor-acceptor molecules, whereby a short side group 
is placed some distance from the donor or acceptor end of the main chain, and a 
linker of about equal length is placed between that side group and another located 
near the other end of the chain. With this configuration, the shape of the eigenstates 
are distorted as displayed in the upper right panel of Figure |o[ generating a large 
change in dipole moment while maintaining a large overlap among the ground state 
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Fig. 24. Molecular form of structures under study. 


and the lowest eigenstates. 



# added single-triple bond inserts, link 

Fig. 25. Variation of f5 X xx and jxxxx for a donor-acceptor structure with a long end group and 
a varying link edge. 


It is now clear that the graph model suggests that a second design rule for 
nonlinear optical structures is to place a defect or side group on a wire somewhere 
between an end and the middle, but not at the end of the chain, as is common 
practice for creating donor-acceptor molecules. Moreover, the presence of potential- 
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altering groups along a chain will do nothing to enhance the optical nonlinear it ies, 
but it can convert a superscaling spectrum into a non-optimum scaling spectrum 
and reducing the nonlinear it ies by many orders of magnitude. Conventional molec¬ 
ular design rules may require some rethinking, but the exact details of how the 
design rules should be stated awaits a many-electron calculation. 


5. Dressed graphs with side groups 

The presence of side groups and defects on a wire, at the same time is important 
to investigate, as it represents a realistic outcome for the design of a quantum wire 
or quasi-linear molecule. Figure ^(i| shows a bare wire and three variants, each with 
both S functions and a prong attached. 





Fig. 26. The bare wire graph and three variants combining prongs and S potentials. 


We examine each in turn. 


5.1. One 5 at one prong 

The upper right panel of Figure p6| shows a one prong wire to which a finite potential 
of the form V(x) = ( g/L)S(x — a) has been placed at the central vertex of a star. 
The solution to this graph has been published Lj. The edge functions match at the 
central ve rtex hi i t. the presence of the delta function changes the flux conservation 
condition l3l§IZ 3 from its canonical form L3 l 2I to 
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F s tar(a, b, c) + (2 g/kL) sin ka sin kb sin kc = 0. (5.1) 


The solutions to Eq. ( [hT| ) for p 
and are described in detail in ref. 


gitive energies (real k) may be found numerically 


5.2. One 5 on one prong 

The lower left panel of Figure ^6] shows a one prong wire to which a finite potential 
of the form V ( x ) = (g/L)S(x — a) has been placed on one prong. This is topologically 
a star graph with a S on one arm. For the 3-star with edges £ 1 , 22 , 23 , the secular 
function F star (zi, 22 , 23 ) for the graph is written in a form that is easy to use with 
the S motif to make bigger graphs. This is 

F s tar(zu Z2, z$) = sin kz\ sin kz2 cos kz% 

+ sin kz\ cos kz 2 sin kzs (5.2) 

+ cos kz\ sin kz2 sin kz% 

For a single unterminated end A, conservation of flux at the central vertex Z pro¬ 
duces the following secular equation relating the amplitudes at the end A and the 
central amplitude Z: 


ZF st ar(z 1 ,Z2,z 3 ) = A sin kz 2 sin kz 3 


(5.3) 


The right-hand side is the net flux through its unterminated vertex A required to 
conserve flux at the central vertex. For the S motif, with one end attached to the 
star at vertex Z and the other end terminated, we haveS 


Fs(ffi;a i,zi) = ^ 




(5.4) 


where we have normalized the £ potential to the total length L = a\ + 21 + 22 + £3 
of the entire graph, as previously discussed. 

We now combine the star motif with the S motif to get the pair of equations for 
the amplitudes A and Z: 


AFs 1 (gi m ,ai,zi) = Zs'mkax (5.5) 

ZF star (z 1 ,z 2 ,z 3 ) = A sin kz2 sin kzs 

Setting the determinant of Eqs. ( ^5|) to zero, we get the secular function for this 
composite graph: 


Fstar— 5 \ — F s i ar (21,22, £3) F § 1 (^1, Ul2i) 
— sin ka\ sin kz2 sin kz% 


(5.6) 








March 18 , 2015 0:21 WSPC/INSTRUCTION FILE flux'interruption 


26 Rick Lytel, Sean M. Mossman, Mark Kuzyk 


We will rewrite this in a form that allows us to extract a common factor sin kz\ 
and remove it from the characteristic equation. This factor is spurious, as it dis¬ 
connects the graph. Using the expressions for the star motif, we get 

Fstar-S ! = FstariZi, z 2 , z 3 )Sp— smkaismkz! 

kL tot 

+ F star (zi, Z2, zs) sin k(a± + z t ) 

— sin ka\ sin kz 2 sin kz% (5.7) 

where L tot = a\ + z\ + Z2 + £3 is the total length of the graph and is unchanged, 
regardless of the position of the S function on the one prong. 

Using Eq. (|ih 2 |), the second and third lines in Eq. ([T?]) may be combined as 

F star (zi,Z2,z 3 )smk(a 1 + 2:1) (5.8) 

— sin ka\ sin kZ2 sin kz% 

= sin kz\ sin k(z2 + Z3) sin k{a\ + z\) 

+ cos kz\ sin kZ2 sin kZ3 sin k(a\ + z\) 

— sin ka\ sin kZ2 sin kz 3 

= sin kz\ sin k(a\ + z\) sin kz2 cos kz^ 

+ sin kz\ sin k(a\ + z{) cos kz<i sin kz% 

+ sin kz2 sin kz 3 [cos kz\ sin k(a\ + z\) — sin tei] 

= sin kz\ sin k(a\ + z\) sin kz2 cos kz% 

+ sin kz\ sin k{a\ + z{) cos kZ2 sin kz% 

+ sin kz2 sin kz^ (cos 2 kz\ — 1) sin ka\ + sin kz\ cos tei cos kz\ 

— sin kz\ sin k{a\ + z{) sin kz2 cos kz% 

+ sin kz\ sin k(a\ + z\) cos kz2 sin kz% 

+ sin kz\ sin kZ2 sin kz% cos k{a\ + z\) 

= sin kz\F star (a\ + z 1 ,z 2 ,z 3 ) 

Combining the last two equations, we get a nice compact form for the secular 
function for the star graph with a S function in one arm: 


F8tar-s 1 (gi m ,ai,zi,Z2,Z3) = sin kzi 


x 


F'star 5 ^2 5 ^ 3 ) 


2gi 

kLtot 


Sin tel + F st ar{cLi +Zi,Z 2 , Z 3 ) 


which is how it must be, because when g — 0, the S function is absent, and the 
graph becomes a 3 -star but with the first prong of length a\ + Z\. 

The nonlinearities of this graph, calculated from a Monte Carlo computation, 
have the extreme values shown in Figure 
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5.3. Two 6 s, one at a prong 

The lower right panel of Figure [26] shows a star with a 5 at its center and one on 
its arm. 

As we’ve seen, for the 3-star with prongs (a, 6 , c) dressed at the center with a S 
function of strength g c /F with L = a + b + c, the secular function becomes l 2 I 

FstarSc = F star (a,b,c ) + (2g c /kL) sin ka sin kb sin kc. (5.9) 


The motif method immediately yield the secular function for the star with a S at 
the center and one on one prong by using F star s c instead of F star in the amplitude 
equations in Eq. (5.5). The result is 


Fstar—8 C 8\ = F star s c (g c , zi,z 2 , z 3 )F Sl (gi m ,ai 3 zi) 

— sin hct\ sin kz 2 sin kz 3 (5.10) 

Finally, we may write this in a simplified form by using Eq. (|5.9| . Factoring out 
a common factor sinfe^i, we get 


Fatar-5 c 5 1 (9i,9c,z 1 ,z 2 ,Z3,a 1 ) = smkz 1 x G; 


where 


+ 

+ 

+ 

where L to t = a\ + z\ + z 2 


Fr s tar—8 c 8i 9ci Z \, Z 2 , £ 3 , U-l) 

Fstar (Zl, *2, Z 3 )-fP~ sin fcoi 
kL tot 

Fstar {CL! +Z 1 ,Z 2 ,Z 3 ) 

c sin k{a\ -\- z{) sin kz 2 sin kz 3 

kL t ot 

sin ka\ sin kz\ sin kz 2 sin kz 3 
kL 2 

+ £3 is the total length of the graph. 


(5.11) 


6. Conclusions 


All of the variants of the bare wire graph shown in Figures ^ [l3|, and 26 have favor¬ 
able geometries (configurations of the edges) to produce very large first and second 


hyperpolarizabilities, as shown in Figure |27[ All have topologies that generate ran¬ 
dom spectra but with fixed gaps set by a quadratic scaling law. For each class of 
graphs, there exist ideal geometries such that the shape of the eigenfunctions may 
be tailored to provide a giant enhancement in their nonlinear optical responses. 

Structures with large intrinsic second and third hyperpolarizabilities necessarily 
have superscaling spectra whose energy E n scales with mode number n as E ~ ttT, 
with fi > 2. This result arises solely from quantum physics and holds for any system 
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Fig. 27. Maximum | /3 XXX \, and maximum and minimum jxxxx for the wire graphs in Figs [ 2 ], [p| 
and|| 


whose nonlinearities may be described by a sum over states. This follows from exten¬ 
sive Monte Carlo analysis of a very large space of spectra and moments constrained 
only by the Thomas-Reich-Kuhn sum rules, which themselves are consequences of 
the validity of the completeness of the eigenfunctions, ie, of the self-adjointness of 
the underlying Hamiltonian, and the validity of the canonical commutation rela¬ 
tions. It is thus a necessary and fundamental design rule that highly active nonlinear 
optical molecules must have superscaling spectra. In a nutshell, this conclusion ex¬ 
plains at once the giant gap between the fundamental limits and all of the measured 
optical nonlinearities to date. But superscaling spectra are not sufficient for large 
optical nonlinearities. 

In this paper, we have shown by use of quantum graph models that systems 
designed with electron flow pathways away from their main chains will create suf¬ 
ficient spatial separation among the ground state and first few excited states to 
generate large differences X 22 and x\\ in their dipole moments required for the 
nonlinear optical response to approach the fundamental limits, provided that the 
geometry of the structure is optimized. In fact, quantum loop graphs have super- 
scaling spectra but very small nonlinearities lZI precisely because their geometry is 
almost maximally non-optimum. But linear and quasi-linear graphs, resembling the 
geometries of conjugated chains with side groups, have near-optimum geometries, 
and when the flux along the main chains is interrupted (ie, discontinuous across 
the chain), the ground and first few excited states are spatially separated, creating 
a large shift in the flow of charge upon a virtual transition to an excited state from 
the ground state, and generating a large nonlinear optical response. 

The two design rules (1) superscaling spectra and (2) spatially separated ground 
and excited states that still maintain good spatial overlap, are heuristics for design¬ 
ing new molecules and nanowires with large hyperpolarizabilities. They do not ex¬ 
plain how to design such molecules, but they do explain why most prior designs fail 
to achieve large values. The heuristics also point toward starting points for design- 
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ing new molecules. For example, structures containing many electrons filling single 
particle states to a Fermi level, will have large responses if the transition moments 
of the electron at the top of the Fermi level to unfilled states above it dominate 
all other transitions from electrons below the Fermi level to states above it. In this 
case, the system behaves as a quasi-one-electron structure, and the results in this 
paper correspond to physically viable models. Such systems, designed as networks 
of structures emulating a quantum graph with a quasi-linear geometry and a side 
prong, will have both superscaling spectra and the transition moments required for 
large hyperpolarizabilities. 
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Appendix A. Off resonance square wells and steps 


The graphs in the main body of the text were wires with zero potential (bare 
graphs) or with S function potentials (dressed graphs). Since the 5 function operates 
at a single point along the wire, its spectrum and that of bare graphs scale with 
energy eigenfunction number in similar ways. For small, positive S strength g the 
two topologies are nearly identical. The effect in each case is to produce a spatial 
shift in the ground and first few excited states whose magnitude depends either on 
the length of the prong or the strength g of the 5 function. This appendix shows 
how wires dressed with potentials of finite spatial extent can produce the same 
spatial shifts required to generate large nonlinear it ies. Dressed graphs with a finite 
potential V(x) are wells or steps. We consider here the square (V(x) = constant ) 
and slant (V(x) ~ x well potentials. The analysis shows behavior similar to those 
of the bare and S dressed wires, with hyperpolarizabilities capable of approaching 
the fundamental limits. This is not surprising, since the dominant feature creating 
the required spatial separation among the lowest eigenfunctions is anything that 
causes a spatial kink in the shape of the eigenfunction at an appropriate location 


on the wire. Figure Al, left, shows the motif graph for the square well. Flux enters 
and exits at the ends and is conserved at the center vertex Z. The edge states are 
written in canonical form that matches flux at the vertices: 


<f>L(x) = 

4>r(x) = 

Conservation of flux yields 


A sin k\ (a — x) + Z sin k\x 
sin k\a 

B sin &2 (x — a) + Z sin &2 {L — x) 


(A.1) 


sm.k 2 (L — cl) 


ZF sa (k 1 ,k 2 -,a,b) = Aki 
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with a secular function 

F sq (k i, &25 b) = ki cos kia sin ^5 + k 2 cos &25sin k\a (A.3) 

and we have set b = L — a. The wave numbers are ki = yj 2 m(E — Vi)/h 2 and 
Vi = (ti 2 /m)g i /L 2 for % = 1,2. 



Fig. Al. Motif for the step on a wire. 


Consider the square well in Fig |A2| , left. The potential for 0 < x < a is 
V{pc) = Vo = (h 2 /m)g/L 2 and zero elsewhere. Note that positive g corresponds 
to a barrier,and all energies are positive. Negative g corresponds to a well, where 
negative energies are possible. Define the wavenumbers for the nonzero potential 
region as K 2 = {k 2 — <a 2 ), where a = ^2g/L 2 . The edge functions are 


= AsiniCr, 0 < x < a (A. 4) 

c/)r(x) = B sin k(L — x), a < x < L 


The secular equation is found by matching the wavefunction and its derivative at 



b 



Fig. A2. Graph for a square well on a wire (left) and a square step on a wire (right). 


the potential boundary. The result is 

k cos k(L — a) sin Ka + K sin k(L — a) cos Ka = 0 


(A.5) 
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It may be seen that this is identical to that obtained by using Eqs. ( |A.2|) and Eq. 
(A.3), with A = B = 0. This may be solved numerically to yield the spectra. The 
edge functions are then used to compute the transition moments. The results are 
displayed in Figures A3 for /3 XXX and A4 jxxx. They are strikingly similar to those 
obtained for a S wire with g < 0, as advertised. 



Position of right edge of potential well 
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Fig. A3. Variation of the first hyperpolarizability Pxxx with the position of the right edge of the 
square step (g > 0) or well (g < 0) potential for twenty values of the potential strength g shown 
in the legend. 


For completeness, we consider the step well in Figure |A2|, ri ght. The flux con¬ 
servation equations are easily found using the motif in Eq. (|A.3| ). They are 

ZiF sq (k\,k 2 m ,a,b) = Z 2 k 2 smk 1 a (A.6) 

Z 2 F sq (ki,k 2 ]c,b) = Z\k 2 sinAqc 

where the wavenumbers on either side of the well are k\ = y/2mE/h 2 and the 
center wavenumber is k 2 = ^/2 m(E — Vo)/h 2 . The secular equation becomes 

F 8 te P (ki,k2 m ,a,c,b) (A.7) 

= F aq (k 1 ,k 2 ; a, b)F 8 q (k 1 ,k 2 ] c, b ) 

— k% sin fciasin k\c — 0 

After a little algebra, a factor sin k 2 b may be removed. Writing F step = sin k 2 b x 
Gstep , we get 

Gstep = \k\ cos k\a cos k\c — k% sin k\a sin k\c\ sin k 2 b 

+ k\k 2 sin (a + c)k\ cos k 2 b (A.8) 
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Position of right edge of potential well 
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Fig. A4. Variation of the first hyperpolarizability f} X xx with the position of the right edge of the 
square step (g > 0) or well (g < 0) potential, again for twenty values of g shown in the legend. 

The spectra may be obtained by solving G step = 0 for the energies. The results are 
similar to those for the square well and are not presented here. 


Appendix B. Off-resonance slant well 

Consider a slant well with infinite barriers at x = 0 and x = L, as shown in Figure 
go is dimensionless. The wave equation takes the form 


Bl. The potential energy (in units where h = m = e = 1 is Vo = (#o/2T 3 )x, where 



0 L 


Fig. Bl. Slant well potential. 


— ^d 2 ip(x)/dx 2 + | V(x) — E] x ) = 0 


(B.l) 
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Define the energy as E = k 2 /2 and a shifted variable £ as 



\ 2 A 

x - L 

\L 3 J 

9o 


(B.2) 


Then the wave equation becomes the Airy equation £^(£) = 0. The solutions 

are 


V’O) = aAi(£) + &£«(£), 


(B.3) 


where £(x) is given in Eq. ( [B.2D and a and b are constants to be determined by the 
boundary conditions that the wave function vanish at x = 0 and x = L: 


where 


aAi(£ 0 ) + bBi(£ 0 ) = 0 
aAii^i) + bBi^i) = 0 

& = -( kLf/gl' 3 
a = 2o /3 [! - (^) 2 /5 o] 


(B.4) 


(B.5) 


From Eq. (R4), we get the secular equation that determines the eigenvalues k n : 

Ai (£o) Bi (£l) - Az(£ L )5i(£ 0 ) = 0, (B.6) 


The solutions to Eq. 


determine the eigenvalues of the slant well. 


Fig B2 plots the hyperpolarizabilities against the potential strength g 0 , as well 
as the ratio of the ground state energy to the potential energy at go. Both the first 
and second hyperpolarizabilities saturate as the well gets higher. 

The eigenfunctions for graphs with go = 0, go = 50, and go = 200 are shown 
below. The effect of a slant well potential is to break the centrosymmetry of the 
go = 0 graph, localizing the ground state of the electron to the left of the well below 
the barrier energy, while pushing the excited states to the right, creating a large 
transition moment and leading to large responses. 

The problem is easily extended to negative potentials, go < 0. The shifted 
position variable becomes 

W) V3r 




[kL) 2 ' 
x+ \ - L 


M 


(B.7) 


and the endpoint values become 


£o = —c kL) 2 /\g 0 \ 2 / 3 
a = -M 1/3 [1 + (kL) 2 /\g 0 \] 


(B.8) 


The secular equation remains unchanged, except for the replacement of the 


endpoint values in Eq. (BJ3) with those in Eq. (IT8). There are negative energy 
solutions for sufficiently large go that are easily discovered by solving the secular 
equation, with k ik. 
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0 50 100 150 200 250 

Potential strength g 0 

Fig. B2. Dependence of /3 X xx and jxxxx on the potential strength go of a slant potential well. 
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Fig. B3. First seven eigenfunctions of the slant well at zero potential, ie, a bare wire. The graph 
is centrosymmetric with zero (3 XX x and the well known value j xxxx = —0.126. 
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